AD-A169  553  METHODE  DES  RAVONS  AVEC  CflLCUL  D' INTEN5ITE  APPLIOUEE  R  1/i, 
LR  PROPAGATION  RNO  <U)  DEFENCE  RESEARCH  ESTA6LISHHENT 
VALCARTIER  (QUEBEC)  C  DION  JJN  06  DREV-4487/86 

F/G  28/14  NL 


UNCLASSIFIED 


UNCLASSII  IED 
U  ML  I  MITE  D  DISTRIBUTiOM 


-vm 

:  r 


CRDV  RAPI’ORi  'MO/  00 
DOSSIER  3MOD  u00 
JUIN  1986 


Ml  I  HODL  DIMIMONS  WH  (  A  I  <  I  I  I  >T  N  I  I  \  Si  |  | 
\  1*1*1  lot  II  \  I  \  I’ KOI*  \(.  \  I  ION  WORM  \l  I 


I)  Dion 


o_ 

o 

o 

UU 
— J 

CO 

l—~ 


•<<  ••  ANI  I  :>l  a  '.OPME  NT  ftHANCH 

.[.c  41. ,,  i^aTiON-H  DEfEN(,f 
A  * A  '  1 A 


-  -4‘  t  ?  OSvHL^PS'tMf-NT 

\s  -(Hi  ~.r  y-i  tn-it  NAT’'.HAlt 
'  ANA'_  A 


i 


Defence  Research  Establishment 
Centre  de  recherches  pour  la  Defense, 

Valcartier,  Quebec 


Canada 


SAMS  CLASSIFICATION 
DISTRIBUTIOM  ILLIMIIEE 


86 


7 


1  4 


1  0  4 


DREV  R-4A07/86 
FILE:  3632D-005 


UNCLASSIFIED 


CRDV  R-4407 / 8 6 
DOSSIER:  3632D-005 


METHODE  DES  RAYONS  AVEC  CALCUL  D’INTENSITE 
APPLIQUEE  A  LA  PROPAGATION  "ANORMALE" 
par 

D.  Dion 


DEFENCE  RESEARCH  ESTABLISHMENT 
CENTRE  DE  RECHERCHES  POUR  LA  DEFENSE 
VALCARTI ER 
Tel:  (418)  844-4271 


Quebec,  Canada 


June/juin  1986 


SANS  CLASSIFICATION 


SAN’S  CLASSIFICATION 
i 


ABSTRACT 


A  ray-tracing  theory  incorporating  the  ray  intensity  calculation 
has  recently  been  developed  at  DREV.  In  this  report,  we  study  the 
applicability  of  this  new  ray  theory  to  radar  propagation.  Using  this 
theory,  we  tried  to  characterize  the  "anomalous"  behavior  of  the  propa¬ 
gation  due  to  the  refractivity  of  the  troposphere. 

This  report  summarizes  the  ray  theory  and  the  validity  of  this 
approach  for  the  study  of  radar  propagation.  A  computer  model  has  been 
developed  to  calculate  ray  trajectories  and  intensities.  Considering 
the  case  of  a  simple  surface  duct,  the  results  have  shown  that  the 
intensity  computation  is  an  effective  and  powerful  way  of  detecting  and 
locating  in  space  such  singular  propagation  phenomena  as  caustics  and 
radio  holes. 


RESUME 


Une  theorie  des  rayons  permettant  le  calcul  d'une  intensite 
associee  a  un  rayon  a  ete  recemment  developpee  au  CRDV.  Dans  ce  rap¬ 
port,  on  etudie  1 'applicabi li te  de  cette  nouvelle  theorie  des  rayons  a 
la  propagation  radar.  A  l'aide  de  la  theorie,  on  a  tente  de  caracteri- 
ser  les  phenomenes  anormaux  de  propagation  causes  par  la  refractivite 
dans  la  troposphere. 

Ce  rapport  rappelle  la  theorie  et  la  validite  de  cette  approche 
pour  l'etude  de  la  propagation  radar.  Un  module  informatique  a  ete 
developpe  pour  calculer  la  trajectoire  de  rayons  et  leur  intensite. 
Considerant  le  cas  d'un  conduit  de  surface  simple,  les  resultats 
montrent  que  le  calcul  de  l'intensite  est  un  moyen  efficace  de  detecter 
et  de  situer  dans  l'espace  les  phenomenes  singuliers  tels  les 
caustiques  et  les  trous  radio. 
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1.0  INTRODUCTION 


La  refractivite  de  l'air  joue  un  role  considerable  dans  la 
detection  electromagnetique;  elle  peut  meme  devenir  dans  certains  cas 
le  facteur  dominant.  Les  references  1  et  2  resuroent  les  differents 
phenomenes  Electromagnetiques  que  provoque  la  variation  de  l'indice  de 
refraction  dans  la  troposphere.  La  reference  2  precise  de  plus  de 
quelle  faqon  ils  sont  relies  aux  conditions  atmospher iques .  Les  pheno¬ 
menes  les  plus  importants  a  signaler  sont  les  conduits  et  les  "trous 
radio".  En  effet,  on  peut  assister  a  une  forte  augmentation  de  la 
couverture  radar  dans  une  certaine  zone  pendant  qu’une  zone  voisine 
presente  un  af faiblissement  marque  de  celle-ci.  De  plus,  un  conduit 
est  constitue  de  rebondissements  multiples  de  la  radiation  qui  s’entre- 
lacent  pour  former  des  caustiques.  II  devient  alors  tres  complexe  de 
tenter  d'evaluer  la  performance  des  systemes  de  communication  ou  de 
detection  dans  ces  conditions.  Ces  phenomenes  etant  tres  frequents  en 
mer,  il  est  d'interet  pour  la  marine  de  rechercher  des  methodes  simples 
perraettant  de  les  caractEriser . 

Des  equations  d ' integration  de  rayons  sur  la  phase  et  une 
methode  de  calcul  de  1 '  "int  insite"  sur  un  rayon  ont  ete  developpees  au 
CRDV  (ref.  3-4).  Cette  theorie  recente  des  rayons  permet  de  repre¬ 
senter  elegamment  le  rayonnement  Electromagnetique  en  milieu  heterogene 
(i.e.  a  indice  de  refraction  variable).  En  particulier,  la  disponibi- 
lite  de  la  phase  et  de  l'intensite  permet  de  sommer  directement  les 
contributions  de  plusieurs  rayons  arrivant  en  un  meme  point.  Par  defi¬ 
nition,  l'intensite  est  la  densite  de  puissance  par  unite  de  surface  en 
un  point  de  l'espace.  Dans  la  theorie  classique  de  1 'electromagne- 
tisrae,  elle  correspond  a  la  grandeur  du  vecteur  de  Pointing.  Dans 
la  theorie  des  rayons,  elle  constitue  une  mesure  directe  de  la  con¬ 
centration  de  rayons  autour  du  rayon  considere.  Jusqu'ici  dans  les 
ouvrages  de  modelisation  de  la  propagation  employant  la  theorie  des 
rayons,  des  approches  indirectes  ou  empiriques  ont  ete  utilisees  pour 
evaluer  l'intensite  du  champ.  Comrae  l'intensite  est  reliee  a  la  den- 
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site  de  rayons  autour  d'un  point,  une  technique  de  calcul  repandue 
consiste  a  emettre  des  rayons  (en  l'occurrence  quatre)  avec  un  leger 
deplacement  dans  les  axes  de  tir  pour  evaluer  1 ' eparpillement  resultant 
des  rayons  sur  une  section  normale  autour  du  point  (ref.  5).  Ceci 
implique  qu'on  doit  repeter  le  processus  en  chaque  point  d'interet.  Une 
autre  solution  est  simplement  de  considerer,  <_n  premiere  approximation, 
un  eparpillement  cylindrique  (en  1/r)  a  l'interieur  d'un  conduit  et 
spherique  (en  1/r-)  a  l'exterieur  (ref.  6-7).  D'autre  part,  on  a  cher- 
che  depuis  les  deux  dernieres  decennies  des  moyens  d'Stendre  la  theorie 
des  rayons  pour  calculer  l'intensite  au  dela  de  surfaces  reflechis- 
santes  et  refractantes  (r£f.  8  a  11).  Les  approches  developpees  sont 
relativement  complexes  et  elles  ne  s'appliquent  pour  la  plupart  qu'a 
des  cas  particuliers . 

Le  but  de  la  presente  £tude  est  de  demontrer  1 'applicabili te  de 
la  theorie,  principalement  en  ce  qui  concerne  le  calcul  de  l'intensite, 
dans  les  conditions  de  propagation  dites  "anormales"  decrites  plus 
haut.  Nous  verrons  qu'il  est  possible  a  l'aide  de  procedes  simples  de 
detecter  et  localiser  les  phenomenes  singuliers  de  la  propagation 
permettant  ainsi  de  caracter iser ,  voire  parametriser ,  le  patron  de 
rayonnement.  Comme  dans  la  plupart  des  cas  les  phenomenes  naissent  a 
la  suite  de  reflexions,  une  methode  permettant  de  poursuivre  le  calcul 
de  l'intensite  a  la  suite  d'une  reflexion  a  du  etre  imaginee. 

Au  chapitre  2,  in  rappelle  brievement  la  theorie  des  rayons 
developpee  au  CRDV  et  au  chapitre  3,  on  discute  de  la  capacite  de  cette 
theorie  de  representer  la  propagation  radar.  Le  chapitre  4  presente 
les  developpements  qu' implique  1 ' applica t ion  de  notre  theorie  3  la 
propagation  t roposphSrique .  La  methode  developpee  qui  permet  de  pour¬ 
suivre  le  calcul  d'intensite  apres  reflexions  y  est  exposee.  Au  chapi¬ 
tre  5,  on  decrit  un  module  informatique  realisant  1 'integration  suivant 
la  theorie  developpee.  Ce  module  poss^de  des  caracteristiques  nou- 
velles  si  on  le  compare  aux  modeles  proposes  dans  la  litterature 
(ref.  5,  6,  12).  II  ne  presente,  entre  autres,  aucune  restriction  sur 
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le  profil  de  l'lndice  de  refraction  n(r,h).  Aussi,  une  terre  ronde  est 
consideree  (quoiqu'il  soit  possible  de  definir  d'autres  formes  de 
terrains)  et  le  calcul  de  la  reflexion  tient  compte  intrinsequement  de 
la  divergence  du  rayonnement  qui  en  resulte.  Enfin,  on  demontre  au 
chapitre  6  la  validite  et  l'efficacite  de  l'approche  en  traitant  le  cas 
d'un  conduit  de  surface  produit  par  un  type  de  profil  de  refractivite 
sinplifi£  mais  d'usage  r£pandu.  Des  techniques  perroettant  d'identifier 
des  regions  singulieres  sont  egalement  presentees  dans  ce  chapitre. 

Ce  travail  a  £te  effectue  au  CRDV  entre  juin  1984  et  aout  1985 
dans  le  cadre  du  NCP  32D05,  Complementarite  des  systemes  optoelectroni- 
ques  et  radars  pour  la  defense  des  navires. 

2.0  THEORIE  DES  RAYONS  ET  DU  CALCUL  DE  L' INTENSITY 


2.1  Equations  d'un  rayon 

En  developpant  le  principe  de  Fermat  (ref.  13),  nous  obtenons 
l'equation  suivante  dite  "eikonal"  qui  constitue  le  fondement  de 
l'optique  geometrique: 


(74>(x,  y,  z))2  =  n2(x,  y,  z)  t2*1] 

ou  n  est  l'indice  de  refraction  carac t£risant  le  milieu  de  propagation 
et  <t>  une  fonction  scalaire  de  position.  Les  rayons  sont  definis  corame 
€tant  les  trajectoires  tirees  or thogonalement  aux  surfaces  4>  ”  cte. 

Ces  surfaces  correspondent  ainsi  3  des  "fronts  d'onde"  optiques. 


Soit  r  un  vecteur  position  dans  l'espace,  on  peut  ecrire: 


SANS  CLASSIFICATION 
4 


ou  Pg  esc  un  vecteur  tangent  au  parcours  d'un  rayon.  En  demontrant 
que : 


ds  d$ 


Blanchard  et  Otis  (ref.  3)  en  d£rivent  une  Equation  dif flrentielle  du 
second  ordre  qui  decrit  le  chemin  "optique”  emprunte  par  le  rayon. 
Elle  s'£crit: 


i>1_  (n2~)  =  Vn  [2.4 

d$  dip 

Cette  expression  a  la  particularity  d'utiliser  p  comme  variable  d'inte- 
gration  au  lieu  de  s:  le  parcours  geometrique.  Elle  peut  etre  trans¬ 
form^  en  deux  Equations  du  premier  ordre  pour  permettre  son  integra¬ 
tion  sur  ordinateur.  On  aura  alors: 


-► 

dr 

d4> 


2.3a 


— —  -  (Vn  -  2n2(Vn*P)p)/n3  [2.3b] 

d  y 

2.2  Equations  de  l'intensite 


En  poursuivant  l'approche  de  Ugincius  (r£f.  14),  Blanchard 
(ref.  4)  propose  une  methode  pour  evaluer  la  variation  de  l'intensite 
causae  par  1 ' €parpi llement  geometrique  des  rayons  en  milieu  het£rogene, 
au  cours  de  1 ' Integration  des  rayons. 

Ugincius  difinit  deux  nouveaux  vecteurs,  a  et  b,  qui  lors  de 
1 ' i ntegrat ion  propageront  1 ' inf orma t ion  sur  la  dispersion  des  rayons 
voisins.  Ces  vecteurs  sont: 


a 


d  r 


[2.6] 
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et 


b  = 


dr, 

d?1 


2.7 


1 

j 


ou  u  et  0  sont  des  parametres  quelconques  d ' ident i fication  des  rayons 
dans  l'espace.  Par  exemple,  dans  le  cas  ou  des  rayons  emaneraient  d'un 
plan,  le  couple  (u),0)  pourralt  etre  les  coordonnees  cart£siennes  (x,  y) 
des  points  d' intersection  des  rayons  et  du  plan. 


Ainsi,  a  et  b  detiennent,  S  un  <£  donne  sur  un  rayon  (w^,  6^), 
de  1 ' inf ormat ion  sur  l'^cart  de  position  qui  existerait  dans  les  deux 

+  -f 

directions  a  et  b  entre  un  rayon  lance  a  (u  +  d^,  0  +  dO)  et  la  posi- 
tion  r  ($)  sur  ce  rayon. 


On  montre  dans  la  r£f.  4  que  l'intensite  le  long  d'un  rayon  est 
alors  decrite  par: 


1(0  =  [n(a  x  b)*p!  U  [2.8] 

|  £ 

En  pratique,  lorsque  les  rayons  se  croisent,  a  la  suite  d'une 
reflexion  ou  en  traversant  une  caustique,  a  ou  b  changeront  de  direc¬ 
tion,  done  de  signe.  On  assiste  en  definitive  a  un  renversement  de 
l'onde  sur  elle-meme.  Pour  eviter  d'avoir  un  1(0  negatif,  il  est 
preferable  de  considerer  a  x  b|  au  lieu  de  (a  x  b)  dans  les  calculs. 
Cette  question  de  renversement  d'onde  sera  rediscutSe  au  chapitre 
su i vant . 


Les  vecteurs  a(0  et  b(0  sont  decrits  par  des  equations  diffe- 

rentielles  du  second  ordre  qui  peuvent  etre  represent£es ,  a  la  maniere 

de  [2.5],  sous  la  forme  de  deux  equations  du  premier  ordre  pour  des 

-► 

raisons  de  calcul  nuraerique.  On  a  pour  a: 


da 


=  (a  -  2n(7n* a)p ) /n' 


2.9a] 
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et  —  =  (n(H  (r)*a)  -  (Vn*a)Vn)/n2 
d  y  n 


ou  H  (r)  est  le  Hesslen  de  n  2  r.  II  s'ecrit: 
n 


2.9b] 


Hn(r) 


32n 

32n 

32n 

sT2 

3y3x 

3z3x 

32n 

32n 

32n 

3x3y 

3z3y 

32n 

32n 

32n 

3x3  z 

3yo«: 

3z2 

De  la  mEme  faqon  pour  b,  on  aura: 


—  =  (6  -  2n(Vn*b)p)/n2 
d: 


et  — L-  =  (n(H  (r)*b)  -  (Vn* b) Vn) /n2 . 


[2.1  la] 


2.11b] 


Les  Eq.  2.5,  2.9  et  2.11  forraent  alnsl  un  systeme  de  18  Equa¬ 
tions  dif ferentielles  du  premier  ordre  decrivant  dans  l'espace  l'evolu- 
tion  du  parcours  des  rayons  et  de  leur  intensite  propre  causEe  par  la 
ref ract ivite.  Le  milieu  de  propagation  est  defini  par  une  fonction 
d'indice  de  refraction  n(r)  qui  peut  etre  quelconque. 

3.0  VALIDITE  DE  LA  THEORIE  DES  RAYONS 


Les  Equations  (ou  lois)  de  Maxwell  constituent  la  thEorie  fonda- 
mentale  pour  decrire  la  radiation  d'energie.  Elies  presentent  le 
phenomeme  sous  la  forme  de  champ  et  d'onde  Electro-magnet iques . 
L'optique  geometrique,  pour  sa  part,  ne  se  penche  pas  d'aussi  pres  sur 
la  nature  merae  du  phEnomene.  En  particulier,  elle  fait  abstraction  du 
concept  d'onde  qui  nous  permet  d'expliquer  et  de  decrire  1 ' inter ference 
et  la  diffraction,  et  de  representer  les  effets  de  la  polarisation. 
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Elle  apparatt  done  incomplete  sur  ce  plan.  Par  contre,  vu  sa  simpli¬ 
city  il  est  grandement  avantageux  d'y  faire  appel  le  plus  possible.  II 
importe  done  de  dlterminer  dans  quelle  mesure  l'optique  g^ometrique 
peut  representer  la  propagation  radar. 

Q 

Depuis  le  debut  du  XX  slide,  on  a  cherche  un  lien  entre  la 
theorie  de  Maxwell  et  celle  des  rayons;  soit  pour  en  arriver  a  des 
equations  semblables  a  celles  de  Maxwell  en  developpant  le  principe  de 
Fermat  (ref.  13),  soit  en  derivant  liquation  eikonal  a  l'aide  de  la 
theorie  de  1 'electromagnetisme. 

Suivant  cette  seconde  voie,  Born  et  Wolf  (ref.  15)  ont  adopte 

les  expressions  suivantes  pour  representer  les  champs  electrique  E  et 
-► 

magnetique  H: 


E  =  E  ej(ko$(r)  '  Ut) 
r 


H  =  H  ej(V(r)  "  ut) 

r 


ou  appele  "nombre  d'onde"  est  egal  3  2n /A .  Ces  auteurs  stipulent 
que  ces  expressions  representent  bien  le  champ  produit  dans  le  vide  par 
un  dipole  monochromatique  et  mesure  a  une  distance  de  plusieurs  lon¬ 
gueurs  d'onde  avec  normalisation  du  moment  dipolaire.  Ils  demontrent 
que,  dans  le  cas  general  ou  E^_  est  complexe,  la  polarisation  pouvant 
etre  elliptique,  on  parvient  I  l'equation  eikonal  en  solutionnant  les 
Equations  de  Maxwell  seuleraent  dans  les  conditions  suivantes: 


-  k  >>  7  x  E 

o  r 


-  k  >>  7  x  H 

o  r 


-  k  >>  E  • 7 lnt  +  7* E 


[3.3b 


[3.3c] 
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-  k  >>  H  • V lnp  +  V*H  [3.3d] 

or  r  L  J 

c 'es t-a-dire ,  a  la  condition  que  la  longueur  d'onde  X  8oit  petite 

devant  le  rotationel  de  E  et  de  H  ,  la  variation  des  caracteristiques 

r  r  n 

du  milieu  (e  et  u )  et  la  variation  du  champ.  C'est  pourquoi  il  est 
generalement  admis  que  l'optique  geometrique  n'est  qu'une  approximation 
qui  n'est  valide  que  lorsque  X  *  0  (ou  kQ  -*■  “).  Ainsi,  en  effet,  la 
longueur  d'onde  et  par  consequent  le  concept  d'onde  s'eteint. 

Cependant,  dans  le  cas  particulier  ou  l'on  considere  une  polari- 
sation  lineaire  et  invariante  de  l'onde  (avec  reel),  on  peut  obteni^ 
l'€quation  eikonal  par  un  d£veloppement  semblable  sans  aucune  approxi¬ 
mation  sur  kQ !  Le  developpement  est  presente  3  l'annexe  A.  Est-ce 
done  alors  qu'une  question  de  polarisation  ou  encore  de  permittivite 
(voir  l'annexe)?  Prenons  garde  de  conclure  trop  rapidement.  II 
advient  que  l'expression  tres  repandue  et  choisie  ici  pour  decrire  les 
champs  E  et  H  (eq.  3.1  et  3.2)  ne  represente  en  realite  qu'une  classe 
tres  restreinte  de  champs  €lectromagnetiques;  en  particulier  en  sup- 
posant  que  E^  est  Independant  de  kQ.  Cette  forme  est,  en  fait,  la 
generalisation  de  l'onde  plane  ou  encore  spherique  en  remplaqant  Ef  par 
E^/r  (ou  r  est  la  distance  parcourue).  Elle  permet  de  les  representer 
parf aitement .  On  peut  ainsi  en  toute  legitimite  emprunter  la  theorie 
des  rayons  dans  la  mesure  ou  1 'approximation  d'onde  plane  ou  sherique 
est  valide.  Aussi,  dans  la  pratique  on  applique  tres  souvent  le  modele 
de  l'onde  plane  pour  simplifier  (ou  souvent  meme  pour  nous  permettre  de 
l'aborder  tout  simplement)  le  traitement  d'un  phenomene  electromagne- 
tique  dans  des  situations  ou  l'onde  n'est  plus  nScessairement  purement 
plane.  C'est  selon  cette  approche,  par  exemple,  qu'ont  ete  derivees 
les  expressions  du  coefficient  de  reflectivite  de  la  radiation  contre 
un  dielectrique  utilise  dans  les  calculs  de  propagation  (ref.  24). 
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Nous  jugeons  alors  preferable  de  qualifier  I'onde  de  "quasi- 
planeite”  ou  “quasi-sphericite“  pour  denoter  les  cas  ou  l'emploi  de 
l'optique  geometrique  est  justifiable.  La  quasi-planeite  ou  sphericlte 
s'estime  a  l'echelle  de  I'onde  dont  la  dimension  se  mesure  en  nombre 
d'onde  (kQ)  ou  longueur  d'onde  (X)*  Ces  caracteristiques  de  I'onde 
s'effondrent  lorsque  le  milieu  ou  l'amplitude  du  champ,  ou  les  deux 
varient  de  faqon  appreciable  en  deqa  de  ces  unites.  C'est  le  cas  lors 
de  la  formation  de  caustiques  ou  l'intensite  de  la  radiation  varie 
considerablement  1  l'interieur  d'une  longueur  d'onde.  Bien  qu'il  ne 
soit  pas  deraisonnable  de  continuer  a  decrire  a  l'aide  de  rayons  le 
parcours  de  la  radiation,  il  n’est  plus  possible  de  tenter  d'evaluer 
son  amplitude.  Tomljanovich  (r£f.  16)  estime  que  l'optique  geometrique 
demeure  valide  jusqu'a  environ  10  ou  20  longueurs  d'onde  autour  d'une 
caustique.  Plus  pres,  on  doit  faire  appel  a  la  theorie  classique  de 
1 'electromagnetisme.  Plusieurs  auteurs  ont  etudie  ce  probleme  du 
calcul  de  l'energie  au  voisinage  des  caustiques.  Les  references  17  et 
18  sont  parmi  les  plus  connues.  L'intensite  montrerait  de  rapides 
oscillations  (battements)  en  traversant  une  caustique  due  aux  franges 
d 'interferences  de  plus  en  plus  rapprochees.  Cette  question,  qui 
constitue  en  elle-meme  un  domaine  de  recherche,  deborde  du  cadre  du 
present  rapport. 

II  est  possible  d'etendre  la  theorie  classique  des  rayons  pour 
representer  le  phenomene  d ' i nter f erence  qui  est  1 'augmentation  ou  la 
diminution  de  l'amplitude  resultant  de  la  superposition  de  rayons  due 
aux  trajets  multiples.  Comme  1' integration  s'effectue  directement  en 
relation  avec  $ ,  on  y  parvient  en  concevant  le  chemin  optique  du  rayon 
subdivise  en  unites  d'onde  auxquelles  on  associe  une  forme  sinusoTdale; 
la  phase  de  I'onde  6tant  donn6e  par  Un  rayon  n'est  done  plus 

uniquement  qu'un  fil  unid imensionne 1  mais  il  possede  desormais  une 
orientation  en  plus  de  sa  direction.  Dans  notre  modele,  cette  orienta¬ 
tion  nous  est  donnee  par  les  vecteurs  a  et  £.  L'intensite  de  la  radia¬ 
tion  a  1 'intersection  de  rayons  est  alors  obtenue  par  la  sommation 
coraplexe  (amplitudes  et  phases)  des  n  rayons  impliques.  On  aura: 
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IT*  =  Reel  {Aje^o^l  +  A2ejko^2  +  ...  +  AneJko*n}  [3.4] 

Le  gain  en  amplitude  risultant  variera  entre  0  et  2n;  ce  qui  correspond 
a  un  gain  en  intensite  variant  de  0  3  4n.  Lorsqu'un  rayon  traverse  une 
caustique,  on  s'aperqoit  qu'il  est  tourne  sur  lui-meme.  Ceci  est 
indique  dans  notre  modele  par  le  changement  d ' orientation  des  vecteurs 
a  ou  b.  On  assisterait  done  3  un  dephasage  de  tt/2.  Les  references  2 
et  3  predisent  en  effet  ce  dephasage.  II  est  a  remarquer  que  le  meme 
phenomene  est  observe  lors  d'une  reflexion. 

La  theorie  des  rayons  ne  permet  cependant  pas  de  reprSsenter  la 
diffraction  qui  explique  la  presence  d'energie  dans  les  zones  ombra- 
gees.  Cette  theorie  montre  en  effet  des  transitions  nettes  entre 
presence  et  absence  d'energie.  Par  exemple,  un  trou  radio  est  claire- 
ment  delimite  et  il  definit  une  region  de  l'espace  ou  l'^nergie  est 
totalement  absente.  On  sait  cependant  qu'en  realite  il  sera  impregne 
du  champ  de  diffraction.  Comme  l'intensite  de  ce  champ  diminue  avec  la 
longueur  d'onde,  on  est  amene  3  considerer,  du  moins  sur  ce  point, 
l'optique  geometrique  comme  un  cas  limite  (lorsque  que  X  -►  0)  de  la 
theorie  de  l'electromagnetisme.  Toutefois,  dans  la  plupart  des  cas  en 
pratique,  le  champ  non  nul  mesure  dans  un  trou  radio  provient  principa- 
lement  des  fuites  caus^es  par  les  irregularites  non  considerees  de  la 
surface  r€f l^chissante  et  du  milieu. 

4.0  APPLICATION  DE  LA  THEORIE  A  LA  PROPAGATION  TROPOSPHERIQUE 

Ce  chapitre  presente  les  d^veloppements  qu'entraTne  l'applica- 
tion  de  la  theorie  (decrite  au  chapitre  2)  3  la  propagation  troposphe 
rique.  Nous  determinerons  les  conditions  initiales  requises  pour 
1 ' integration  des  equations  differentielles  et  nous  exposerons  une 
methode  qui  permet  de  poursuivre  1 'integration  apres  reflexions  de 
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Dans  la  troposphere,  nous  considerons  que  le  milieu  de  propaga¬ 
tion  est  homogene  en  azimut.  La  fonction  d'indict  de  refraction 
decrivant  le  milieu  ne  depend  alors  que  de  l'elevation  (h)  et  de  la 
distance  radiale  a  l'antenne  (d).  Dans  un  referentiel  cartesien  dont 
l'origine  est  situ£e  a  l'antenne,  elle  peut  s'ecrire: 

n(r)  =  n(d,h)  ou  d  =  V  x2+  y 2  [4.l] 

et  h  =  z 

Dans  ce  cas,  nous  n'avons  qu'a  determiner  le  rayonnement  uniquement 
dans  un  plan  transversal  au  sol  pour  representer  entierement  la  propa¬ 
gation. 


En  pratique,  l'indice  de  refraction  qui  est  mesure  a  l’aide  de 
ballons  sondes  n'est  en  fait  donne,  la  plupart  du  temps,  qu'en  fonction 
de  l'elevation.  Cette  simplification  est  particulierement  justifiee  en 
raer  ou  generalement  le  gradient  de  l'indice  selon  d  est  environ  100 
fois  inferieur  a  celui  selon  h.  Ceci  explique  pourquoi  dans  1 'ensemble 
les  modeles  de  propagation  existants  ne  considerent  que  la  variation  de 
n  en  elevation.  Aussi,  la  simplification  de  la  fonction  de  refracti- 
vite  reduit  considerablement  la  complexite  et  la  lourdeur  des  calculs 
quelle  que  soit  la  theorie  adoptee  pour  representer  la  propagation. 
L'etude  d'un  profil  plus  detaille  se  revelait  impraticable  il  y  a 
quelques  annees  seulement. 

Le  gradient  selon  d  devient  cependant  significatif  pres  des 
cotes  en  raison  des  influences  topographiques  sur  la  circulation  de 
l'air  a  faible  altitude  (ref.  19).  Notre  modele  de  n  nous  permet  done 
d'etudier  les  effets  encore  peu  explores  de  cette  composante  radiale  du 
gradient . 
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4 . 1  Conditions  Inltlales 

La  valeur  initiale  des  variables  d 'integration  des  Equations 
dif ferentielles  constituant  notre  modele  depend  du  patron  de  radiation 
spatial  defini  a  l'emission.  Nous  considererons  le  cas  d'une  emission 
spherique,  dite  isotropique,  qui  sert  de  reference  dans  le  domaine  des 
antennes.  Le  rSfSrentiel  utilise  est  montre  a  la  fig.  1  et  l'antenne 
est  situee  &  l'origine  du  systeme  de  coordonnees.  Soit  une  sphere 
inf initesimale  de  rayon  ds  centree  a  l'origine,  un  point  r  sur  la 
sphere  est  represents,  en  cart£sien,  en  fonction  des  parametres  d'iden' 
tification  u  et  0  selon: 


FIGURE  1  -  Le  ref£rentiel  utilise 
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+  +  -►  ^ 

On  obtient  alors  p,  a  et  b  a  l'aide  de  leur  definition  respective 
donnee  aux  eq.  2.5,  2.6  et  2.7.  On  aura: 


dr 


1  O  1  1 

=  _  =  -■  (cosu  cosO,  sinu,  cosu  sinO)  [4.3  i 

ds  n,  ds  n  , 
ds  ds 


ou  n,  =  n(r  ), 
ds  ds 


ds 


=  ds(~sinu  cosO,  cosj,  -sinu  sin6) 


et  b 


ds 


=  ds(-cosu  sinO,  0,  cosu  cosO) 


[4.A] 

[4.5  1 


Les  valeurs  des  parametres  a  I’origine  s’obtiennent  en  proc£dant 
au  passage  a  la  limite  ds  -*-0.  On  a  ainsi, 


rQ  =  (0,  0,  0) 


p  -  —  (cosu  cosO,  sinu,  cosu  sinO) 
o  n 


aQ  =  (0,  0,  0) 


bQ  =  (0,  0,  0). 


[4.6a  ] 
[4.6b] 

[4.6c  ] 
[ 4 . 6d  ] 


_  da  d  ,dr.  d  .dr.  d  T 

Coniine  —  ■  —  (-— )  *  —  (•—-)  =  —  P 

dy  d<J  du  du  d$  du 


on  calcule  a  en  l'isolant  dans  l'eq.  2.9a.  On  obtient: 
o 


a  =  n  (-sinu  cosO,  cosu,  -sinu  sinO), 


[4.7 


De  facon  similaire  —  =  — ,  et  en  utilisant  l'eq.  2.11a  on  obtient  pour 
-►  d<f  dO 

e  : 

o 


n  (-sinO  cosu,  0,  cosO  cosu). 


4.8 


+  Q.--I  (DO.  ;r  Or-lxiB—t/cocucr  +  .£> 
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Dans  l'hypothese  d 'homogenei te  azimutale,  les  expressions  pour 

,  a  et  8  peuvent  etre  simplifiees.  En  effet,  en  ne  considerant  qu< 
ooo  ^ 

e  plan  correspondant  a  w  =  0,  on  aura: 


p  =  —  (cosJ,  0,  sinG)  [4. 9a] 

on  1  1 

o 

S  =  (0,  nQ,  0)  [4.9b] 

6q  =  no(-sin0,  0,  cosG)  [4.9c] 

■+ 

II  est  interessant  de  remarquer  que  est  constant  puisque: 

Vn*a  =  H  • a  =  0, 
n 

t  que  l'ecartement  des  rayons  dans  la  direction  de  a  se  resume  alors 
ar  la  simple  expression: 


n 

da  o  r 

df>  n^"(T)" 

. 2  Calcul  de  l'intenslte  apres  reflexion 


4.10] 


Les  Equations  2.5,  2.9  et  2.11  forment  un  systeme  permettant  de 
alculer  la  trajectoire  de  rayons  et  l'intenslte  qui  leur  est  associee 
'application  des  Equations  de  rayons  est  limitee  si  leur  integration 
olt  s'arreter  lorsque  les  rayons  rencontrent  une  discontinuity  du 
ilieu  qui  provoquerait  une  reflexion.  D'autant  plus  que,  comme  nous 
e  verrons  au  chapitre  suivant,  en  general  les  phenom?nes  interessants 
etudier  prennent  naissance  a  la  suite  de  reflexions.  Une  loi  simple 
t  bien  connue  gouverne  la  direction  des  rayons  reflechis,  mais 
u’advient-il  de  notre  calcul  de  l'intenslte?  Comment  les  vecteurs 
,  a  et  8  sont-ils  affectes  par  la  reflexion?  Pour  rendre  plus  com- 


BJ  +• 
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plete  notre  theorie  des  rayons  et  Itendre  son  applicability,  il  est 
necessaire  de  pouvoir  poursuivre  1 ' integrat ion  de  l'intensite  au-delS 
des  reflexions.  Une  methode  numlrique  simple  a  du  ainsi  etre  conque 
pour  permettre  de  prolonger  les  rayons  5  la  suite  d'une  reflexion. 

Les  nouveaux  vecteurs  d'intensite  (a,  b,  a,  3)  dependent  du 
patron  de  rayonneraent  r£fl£chi.  Ils  doivent  etre  calcules  comme  a 
l'origine  a  l'aide  de  leur  definition  respective.  II  faut  done  deter¬ 
miner: 


* 

a 

r 


a 

r 


+  2n  (Vn  *a  )p  , 
r  r  r  r 


[4.11 ' 


4 . 1 2  ] 


b 


r 


dr 


* 


» 


r 


et  6 

r 


=  nr7T 


2n  (Vn*b 
r  r 


>.13] 


>.14 


Pour  demeurer  cohlrent  avec  la  theorie,  on  est  tenu  d'evaluer 

ces  expressions,  comme  il  est  indique,  3  un  meme  4>  (i*e.  en  phase)  afin 

de  representer  ad^quateraent  le  front  d'onde  optique  reflechi.  Dans 

l'hypothese  d ' horaogenei te  azimutale  et  avec  une  surface  ref lechisante 

spherique,  on  salt  a  priori  que  a  et  a  seront  inchanges.  La  diffi- 

>  -► 

culty  est  done  d'evaluer  —  et  —  aux  points  de  reflexion. 

d3  dO 
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Pour  y  arriver  numeriquement ,  on  doit  connaftre  les  points  de  reflexion 
des  rayons  voisins.  Ceci  nous  oblige  3  attendre  que  tous  les  rayons 
soient  r£fl£chis  et  a  proc£der  ensuite  a  une  reemission  des  rayons.  A 
la  reflexion  d'un  rayon  on  recueillera  l'information  necessaire  pour 
caracteriser  le  "rebondissement".  On  appel  "rebondissement"  un  ensem¬ 
ble  de  rayons  reflechis  ayant  une  origine  commune.  L'information 
pertinente  est  pour  chaque  rayon: 

■f 

les  coordonnees  du  point  de  reflexion:  r  , 

-► 

la  direction  du  rayon  une  fois  reflechi:  p^, 
la  longueur  du  chemin  optique  au  point  de  reflexion:  <p  , 
l'intensite  du  rayon  au  point  de  reflexion:  I^_,  et 
l'identificatlon  du  rayon  0O» 

Localeraent,  autour  d'un  rayon  donn£,  on  reconstitue  le  front 
d'onde  reflechi  en  remettant  les  rayons  voisins  en  phase  avec  celui-ci. 
Quoique  plusieurs  rayons  voisins  peuvent  £tre  utilises,  on  peut  a  la 
limite  n'en  utiliser  qu'un  seul  et  alleger  ainsi  considerablement  le 
calcul. 


Le  proc£de  de  remise  en  phase  est  illustre  a  la  fig.  2.  Un 
rayon  refiechi  plus  loin  sur  la  surface  doit  etre  recule  de  dans 

la  direction  oppos^e  a  sa  reflexion  normale  pendant  qu'un  rayon  refle¬ 
chi  avant,  s'il  est  utilise,  doit  etre  avance.  Le  recul  revient  en 
fait  a  reproduire  les  rayons  images  (par  similitude  a  la  reflexion 
contre  un  miroir)  jusqu'au  premier  raccordement  de  phase.  La  propaga¬ 
tion  sous  la  surface  r£ f 14 ch issant e  souleve  un  probleme  de  definition 
de  l'indice  de  refraction.  Le  recul  nous  oblige  3  reconstituer  un 


*  a 
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systeme  referentiel  et  un  profil  de  r£fractivite  inverse,  et  a  transpo- 
ser  a  la  suite  du  calcul  les  resultats  r(<f>r)  et  p(t^)  obtenus  dans  le 
referentiel  original.  En  pratique,  si  un  nombre  limite  de  rayons  rela- 
tivement  rapproches  est  utilise,  le  parcours  du  recul  et  la  profondeur 
de  penetration  deviennent  relativement  petits.  On  peut  alors,  pour 
simplifier  le  calcul  sous  la  surface,  faire  l'hypothese  d'un  gradient 
de  l'indice  constant  et  €gal  a  celui  au  sol. 

-+ 

On  obtient  de  cette  faqon  des  positions  r(y^)  et  des  directions 
(^r)  qui  definissent  un  patron  de  rayonnement  discret  autour  d'un 
rayon  reflechi.  II  est  alors  possible  d'evaluer  numeriquement  les 
derivees  necessaires  enumer£es  plus  haut. 


En  resume,  dans  l'hypothese  d 'homogenei te  azinutale,  les  valeurs 
initiales  des  parametres  a  la  reemission  seront  pour  chaque  rayon: 


r^:  le  point  de  reflexion 


+  + 

P  =  p  -  2(p  *u  )u 
r  i  inn 


4.15 


ou  u  est  un  vecteur  unitalre  normal  3  la  surface  ref lechissante . 
n 


-  ar  -  (0, 


,  0) 


4.16 


-  ar  =  (0,  nr,  0) 


4.17 


dr 


b  =  -j— — 

r  au 


4.18 


dP 


8  =  n  *■  j  p 

r  r  du 


2n  (Vn  *b  )p 
r  r  r  r 


4.19 


ou  n  = 


n(rr). 
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En  terminant  il  est  important  de  souligner  que  le  patron  de 
reemission  se  definit  lors  du  calcul  des  directions  des  rayons  refle- 
chis.  Ces  directions  sont  relives  a  la  g£ometrie  de  la  surface  r£fle 
chissante  et,  comme  elles  sont  calculees  individuellement  et  ponctuel 
lement,  cette  surface  peut  etre  quelconque  voire  meme  Itablie  numeri- 


quement . 


DROITE  D’EQUI-PHASE 


FIGURE  2  -  Reflexion  dos  rayons 


5.0  MODULE  INFORMAT I QUE 


Un  module  informatique  de  base  a  £te  dSveloppe  pour  effectuer 
les  calculs  de  la  trajectoire  et  de  l'intensite  de  rayons  suivant 
l'approche  decrite  dans  les  chapitres  precedents.  Le  r£f£rentiel 
utilise  est  celui  de  la  fig.  1.  L'origine  est  situee  3  l'antenne  et  on 
definit  sous  celle-ci  une  surface  sph£rique  r§f lechissante  repr£sentant 
la  terre.  Les  rayons  seront  emis  selon  0  pour  =  0  en  utilisant  les 
valeurs  initiales  des  parametres  calcul£es  au  chapitre  4. 
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5.1  Les  entrees-sorties 


Les  parametres  d'entr£e  sont  ajustes  a  l'aide  d'une  routine 
interactive  appelee  "RAYIN:EX".  Ce  sont: 

-  la  hauteur  de  l'antenne  par  rapport  S  la  terre  (h  ) 
le  pas  d 'integration  (A<£) 

le  nombre  de  rayons  S  emettre  (Nr) 

l'angle  du  premier  rayon  (0  ) 

min 

-  le  pas  de  l'angle  d'emission  (AG) 

l'arc  terrestre  maximum  sur  lequel  seront  emis  les  rayons 

(S  ) 
max 

-  le  nombre  d ' echant illons  de  position  et  d'intensite  a  conser 
ver  comme  resultats  (Ng). 

Ces  donnees  sont  emmagasinles  dans  un  fichier  (RAYINPUTS)  pour  etre 
lues  par  le  programme  principal. 

Le  module  produit  quatre  fichiers  de  resultats:  HFIL,  WFIL, 
FIFIL  et  RFLFIL.  Les  trois  premiers  renferment  respectivement  les 
echantillons  d'eievation,  d'intensite  et  de  longueur  optique  pour 
chaque  rayon.  Ce  6ont  des  matrices  de  Ng  lignes  et  Nr  colonnes.  Le 
fichier  RFLFIL  contient  pour  sa  part  1 'inf ormation  decrivant  les 
reflexions  rencontrees.  Chaque  ligne  correspond  a  une  reflexion. 

L ' inf ormation  contenue  est  celle  requise  pour  effectuer  la  reflexion 
des  rayons  par  reemission  en  phase  telle  que  decrite  au  chapitre 
precedent . 

5.2  Le  derouleme  des  operations 

Apres  une  phase  d ' ini t ialisa tion  et  de  lecture  des  donnfes,  le 
module  de  propagation  procede  3  1 ' integration  des  rayons.  Celle-ci  est 
r£alis£e  a  l'aide  de  la  m€thode  de  Runge-Kutta  d'ordre  4  (ref.  20).  A 
long  de  leur  course,  les  rayons  franchissent  les  points  d ' £chant i 1 Ion- 
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nage  S^  sur  la  terre  ou  l'on  evalue  leur  elevation  radiale  par 
rapport  a  la  terre,  leur  intensite  et  leur  longueur  optique  G ^ . 

On  procede  par  interpolation  lineaire  sur  le  pas  d 'integration.  Les 
S^  sont  determines  a  partir  des  parametres  d'entree  selon  la  rela¬ 
tion: 


S 


i+1 


S 

max 

~T7~ 

e 


+  S 


.5-1  j 


pour  i  =  1 . N,  -  1  et  S,  =  0. 

e  1 

Les  resultats  produits  sont  accumules  respectivement  dans  les  fichiers 
HFIL,  WFIL  et  FIFIL. 


Pour  chaque  rayon,  1 ' integration  s'arrete  lorsqu'il  atteint  la 
distance  limite  ^  ou  une  reflexion  (lorsqu'il  tente  de  penetrer  a 
l'interieur  de  la  surface  de  la  terre).  Dans  ce  dernier  cas,  on 
caicule  les  parametres  carac ter Isant  le  rebondissement  en  formation 
(voir  chapitre  A),  soit: 


le  point  de  reflexion:  r 

+ 

la  direction  du  ravon  reflechi:  P 

r 

son  intensite  au  point  de  reflexion: 

et  la  longueur  optique  au  point  de  reflexion: 


Ces  valeurs,  jointes  a  des  donnees  d ' ident if i cat  ion  des  rayons  et  du 
rebondissement,  sont  notees  au  fichier  RFLFIL. 


Lorsque  tous  les  rayons  ont  ete  emis,  on  verifie  s'il  y  a  eu  un 

rebondissement.  Dans  1 ' af f i rmat i ve ,  on  procede  5  la  reemission  des 

rayons  reflechis  telle  que  decrite  au  chapitre  A.  Les  deriv£es  num(ri 

ques  servant  a  calculer  les  valeurs  initiales  de  b  et  C  a  la  reemis- 

r  r 

sion  (eq.  A. 13  et  A.1A)  sont  evaluees  a  l'aide  d'un  polynome  de  colla¬ 
tion  d'ordre  2  en  considerant  pour  chaque  noeud  ses  voisins  de  droite 


s  *» 
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et  de  gauche.  Aux  extremites,  on  doit  utiliser  1 'approximation 

lineaire.  Apres  etre  reflechis,  les  rayons,  5  l'exemple  de  ceux  qui 

emanent  de  l'origine,  sont  propages  jusqu'a  ce  qu'ils  rencontrent  la 

distance  maximale  S  ou  un  nouveau  rebondissement •  Le  processus  se 

max 

termine  lorsque  tous  les  rayons  sont  parvenus  a  S 

max 

5 . 3  Structure  du  loglciel 

Le  module  informatique  realise  est  forme  de  quatre  blocs.  Le 
programme  directeur  et  les  routines  principales  sont  contenus  dans 
PROPAK,  UTILPAK  renferme  des  routines  utilitaires  de  calcul,  NRPAK  les 
routines  d'operations  vectorielles  sur  l'indice  de  refraction  et  FNPAK 
les  fonctions  de  refractivite  elles-memes.  Le  contenu  de  chacun  des 
blocs  est  resume  aux  tableaux  I  a  IV.  L'asterisque  accompagnant  le  non 
d'une  routine  signifie  qu'il  s'agit  d'une  fonction  FORTRAN  et  non  d'une 
sous-routine .  Ces  blocs  peuvent  etre  compiles  et  lies  par  la  commande 
JCL  "PR0PEX:T."  qui  produit  le  module  executable  "PROPAG:EX". 


TABLEAU  I 

Bloc  principal  PROPAK 


NOM  DES  ROUTINES 

FONCTION 

PROPAG 

Programme  directeur 

EMITRAY 

Coordonne  l'emission  d'un  rayon 

RAYSTP 

Effectue  un  pas  d 'integration  par  la  methode  de 
Runge-Kutta 

PROPFS 

Appele  par  RAYSTP  pour  evaluer  les  six  equations 
de  la  theorie  des  rayons 

REFLX 

Appele  par  EMITRAY  pour  traiter  une  reflexion 

STOREZW 

Appele  par  EMITRAY  pour  calculer  et  porter  a  des 
flchlers  respectifs  Hj,  Wj  et  v p  si  le  rayon  a 
frauchi  la  distance  d'echantillonnage  courante 

RE EM IT 

Une  fois  tous  les  rayons  emis,  realise  la 
reemission  des  rayons  reflechis 
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TABLEAU  II 

Bloc  utilitaire  UTILPAK 


NOM  DES  ROUTINES 

FONCTION 

*  H 

Calcule  a  r  la  hauteur  Bur  la  normale  a  la  terre 

*  S 

Calcule  S  r  la  distance  terrestre  equivalente 

SCAPROD 

Effectue  le  produit  scalaire  entre  deux  vecteurs 

VECPROD 

Effectue  le  produit  vectoriel  entre  deux  vecteurs 

*  DEGgRAD 

Conversion  de  radians  a  degres 

*  RADgDEG 

Conversion  de  degres  a  radians 

DIRVECT 

Retourne  la  grandeur  et  le  vecteur  unitaire 
correspondant  au  vecteur  d'entr£e 

TABLEAU  III 

Bloc  NRPAK  d'operateurs  vector lels  sur  n 


NOM  DES  ROUTINES 

FONCTION 

RPRIM 

Calcule  la  coordonnle  r',  dans  le  referentiel  ou 
est  decrit  n(d,h),  correspondant  a  la  position 
r  dans  le  referentiel  de  base 

*  FNR 

Donne  la  valeur  de  n  au  pt.  r  en  appelant  RPRIM 
et  FN0 

GRADIENT 

Calcule  Vn  au  pt.  r 

HESSIAN 

Calcule  le  produit  Hn* V  au  pt.  r;  V  £tant  un 

vecteur  quelconque  specific  coiame  argument 
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TABLEAU  IV 


Bloc  des  fonctlons  de  refractivite  FNPAK 


NOM  DES  ROUTINES 

FONCTION 

*  FN0 

Fonction  n(d,h;  represente  dans  son  referentiel 

*  FN1 

0 

Fonction  —  dans  le  referentiel  de  n 

3h 

*  FN2 

a  2 

Fonction  S. 

3h2 

6.0  ETUDE  DU  CAS  D'UN  CONDUIT  DE  SURFACE  SIMPLE 


Dans  le  but  de  verifier  1 'applicabili te  du  modele  theorique 
presente  et  de  demontrer  les  possibilites  du  module  informatique  ele- 
mentaire  qui  en  decoule  directement,  on  a  considere  le  probleme  de  la 
caracterisation  de  la  propagation  dans  des  conditions  de  conduits 
tropospheriques  de  surface. 

Pour  des  fins  de  comparaison  avec  des  modeles  ou  des  etudes 
connues,  un  profil  d'indice  de  refraction  approxime  de  segments  de 
droite  a  ete  choisi.  II  est  montre  a  la  fig.  3.  Le  graphique  pre¬ 
sente,  en  fonction  de  l'elevation  h,  la  variable  N  appelee  indice  de 
refraction  modifie  ou  encore  ref ractivi te;  elle  est  relive  a  l'indice 
de  refraction  par  la  relation: 

N  -  (n-1)  x  10~6.  [6.1] 

Ce  profil  simplifie,  qui  neglige  la  seconde  derivee  et  provoque  une 
discontinuity  de  la  premiere,  est  encore  le  plus  utilise  dans  les 
Etudes  de  propagation  dite  guid£e  (ref.  5,  12  et  21).  Deux  raisons 
justifient  ce  choix.  D'une  part,  les  modeles  de  propagation  developpes 
jusqu'ici,  soit  a  partir  de  la  theorie  des  modes  (ref.  21)  ou  des 
rayons  (ref.  5),  ne  permettent  pas  pour  la  plupart  de  considerer  un 
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type  de  profil  plus  precis.  D'autre  part,  une  theorie  populaire  relie 
la  canalisation  d'energie  a  la  presence  de  "mur"  atmospher ique  cause 
par  la  discontinui te  de  n'. 

Les  grandeurs  des  parametres  du  profil  utilise  dans  notre  etude, 
soit  les  2  pentes,  la  hauteur  du  conduit  et  la  valeur  de  N(0),  sont 
celles  utilisees  par  Boithias  (r£f.  22).  Selon  l'auteur,  ce  sont  des 
valeurs  assez  representatives  dans  des  regions  ou  la  propagation  guidee 
est  persistante. 

6. 1  Trace  des  trajectoires  des  rayons 

Analysons  brievement  les  diagrammes  de  rayonneraent  obtenus  pour 
differentes  hauteurs  d'antenne  (fig.  4,  5  et  6).  La  couverture  est  de 
1.125°  et  les  rayons  sont  emis  a  des  intervalles  de  .075°.  Meme  avec 
un  nombre  restreint  de  rayons  (en  l'occurrence  15),  on  parvient  a 
distinguer  netteraent  les  regions  singulieres:  trous  radio  et  causti- 
ques  se  rejoignant  a  des  cuspides.  Ces  resultats  sont  conformes  ii  ceux 
de  Boithias. 

A  partir  de  ces  resultats,  on  peut  degager  des  observations 
pertinentes  sur  la  formation  de  conduits  de  surface  simples: 

-  les  cuspides  et  le  sommet  des  caustiques  principales  se 
situent  relativement  3  la  hauteur  de  l'antenne, 

la  distance  entre  les  bonds  diminue  avec  la  hauteur  de 
1 'antenne , 

la  distance  entre  l'antenne  et  la  premiere  cuspide  est 
approxima t i vement  egale  3  la  distance  entre  deux  bonds, 

un  trou  radio  prend  origine  a  la  distance  de  la  premiere 
cuspide  et  3  la  hauteur  du  conduit. 
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De  plus,  on  salt  que  l'effet  de  canalisation  d'energie  est  inexistant 
lorsque  l'antenne  est  situee  au-dessus  de  la  hauteur  du  conduit.  Ce 
parametre  qu'est  la  "hauteur  du  conduit”  est  grandement  utilise  dans  la 
caracter isat ion  d'un  conduit.  II  est  bien  defini  dans  cet  exemple-ci 
puisqu'il  coincide  avec  la  hauteur  de  la  discontinue  de  N’.  Cependant, 
il  devient  plus  difficile  3  determiner  dans  le  cas  general  ou  N 
provient  de  mesures  exp£ rimentales  discretes.  Aussi,  on  en  retrouve 
plusieurs  definitions  mathematiques  dans  la  litterature.  Jeske 
(ref.  23)  relie  la  hauteur  du  conduit  a  la  temperature,  a  la  vitesse  du 
vent  et  a  la  temperature  de  la  mer  (ainsi  a  la  stabilite  de  l'atmo- 
sphere);  cette  definition  est  probablement  la  plus  utillsee. 


FIGURE  A  -  Diagramrae  de  rayonnement 


avec  h  = 
a 


50  m 
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Cette  technique  du  trace  des  trajectoires  des  rayons  est  bien 
connue.  Examinons  done  plus  attentivement  1 ' information  nouvelle 
apportee  par  l'addition  du  calcul  de  l'intensite  pour  d£couvrir  dans 
quelle  mesure  elle  peut  contribuer  a  mieux  caracteriser  la  propagation. 
A  priori,  on  devrait  s'attendre  intuitivement  que  cette  mesure  de 
l'intensite  nous  procure  une  grande  facilite  2  deceler  de  faqon  numeri- 
que  la  presence  de  zones  singulieres  puisque  l'intensite  devrait  d'une 
part  former  des  pics  aux  cuspides  et  caustiques,  et  d'autre  part  subir 
un  af faiblissement  marque  aux  trous. 


6.2  Comportement  de  l'intensite  des  ravons  pres  d'un  trou  radio 


Considerons  le  cas  de  la  fig.  6  (h„  =  190m).  Deux  trous  radio 
sont  presents:  le  premier  a  l'origine  et  le  second  a  la  premiere 
cuspide.  Puisque  la  localisation  de  la  cuspide  determine  l'origine  du 
second  trou,  interessons-nous  done  a  celui  forme  a  l'antenne. 


A  diverses  distances  terrestres  S^^  de  l'antenne,  on  examine 
les  valeurs  d'intensite  sur  les  rayons  (celles-ci  6ont  accumul£es  lors 
du  calcul  au  fichier  WFIL  (voir  chapitre  5)).  On  remarque  qu'en 
s'approchant  ou  en  s'^loignant  du  trou,  l'intensite  ne  varie  que  tres 
tres  faiblement:  environ  .0001  dB  par  rayon.  Le  tableau  V  presente 
les  valeurs  d'intensite  des  rayons  pres  du  trou  pour  trois  distances 
d 'echantillonnage.  Cependant  pour  suffisamment  grand,  3  30.4  km 
par  exeraple,  on  observe  un  saut  toujours  faible  en  valeur  absolue,  en 
l'occurrence  .0012  dB  (lequel  correspond  3  une  variation  d'intensite  de 
1.00027),  mais  qui  est  significatif  devant  les  hearts  d'intensite 
existants  entre  les  autres  rayons.  On  obtient  ainsi  un  moyen  d'identi- 
fier  les  rayons  limites  au  trou.  Les  r£sultats  du  tableau  V  peuvent 
6erabler  en  contradiction  avec  le  diagramme  de  rayonnement  qui  montre  un 
espaceraent  progressif  assez  marque  des  rayons.  II  faut  cependant  £tre 
tres  prudent  dans  1 ' interpretat ion  de  la  figure  6.  On  doit  tenir 


SANS  CLASSIFICATION 
29 


compte  de  l'importante  difference  dans  les  £chelles,  I'abcisse  etant 
environ  150  fois  plus  etendue  que  l'ordonnee.  Cette  forme  de  graphlque 
perraet  d'accentuer  les  ph£nomenes  qui  nous  int£ressent  mais  elle  peut 
etre  tronpeuse. 

La  variation  de  l'intensite  causee  par  la  presence  du  trou  S 
1'origine  est  grandement  submergee  par  1 'eparpl llement  naturel  en  1/r2. 
II  pourraft  done  £tre  preferable  dans  ce  cas  de  considerer  plutot  la 
valeur  de  l'intensite  normalisee  a  celle  qu'on  auralt  en  milieu  homo- 
gene:  I.  Elle  est  presentee  a  la  troisieme  colonne  du  tableau  V. 

Celle-ci  montre  la  faible  divergence  par  rapport  au  milieu  homogene.  Vu 
les  grandeurs  en  jeu,  afin  d'exploiter  adequatement  ce  parametre,  on 
doit  tenir  compte  de s  erreurs  inherentes  aux  calculs  num£riques.  On 
est  ainsi  amene  a  determiner  un  seuil  d’erreur  en-de?a  duquel  I  n'est 
plus  signif icatif .  Dans  cet  exemple  il  a  ete  fixe  empi riquemen t  a 
10-u  dB.  Bien  que  l'erreur  (en  fait  son  ecart  type)  crolt  avec  la 
longueur  du  rayon,  on  estime  qu’il  n’excedera  pas  ce  seuil  a  l'inte- 
rieur  des  limites  considerees.  Des  tests  effectues  en  double  precision 
ont  permis  de  corroborer  cette  hypothese.  Comme  le  montre  le  tableau  V 
ce  nouveau  parametre  met  plus  en  evidence  les  variations  (relatives) 
d’intensite.  On  s'aperqolt  que  ce  n'est  pas  l'intensite  (ou  meme 
l'intensite  relative)  elle-meme  qui  est  revelatrice  pour  situer  ce  trou 
mais  plutot  les  differences  d'intensite  ou,  en  d'autres  termes,  sa 
der i vee . 


S  AN S  C  LA  SSI  FI  CAT  I  ON 
30 


TABLEAU  V 


Nesures  d'intenslte  do  rayons  autour  du  trou  radio  a  l'orlgine 
_ (N'S:  Valour  conslderee  non-signi f icat Ive  ) _ 


Rayon 

I(dB) 

Al(dB) 

I(dB) 

6 

-80.9844 

9.16  x  10~->  (NS) 

7 

-80.9845 

.0001 

9.16  x  10~5  (NS) 

S  =  11.2  km 

8 

-80.9846 

.0001 

~0  (NS) 

9 

-80.9847 

.0001 

~0  (NS) 

6 

-86. 3607 

.74  x  1 0~  ' 

7 

-86.3608 

.0001 

.70  x  1  0"  J 

S  =  2o.fi  km 

8 

-86.3613 

.0007 

~0  (NS) 

9 

-86.3616 

.0001 

6.10  x  10"5  (NS) 

6 

-89.6565 

1.1  x  10" 

7 

-89.6566 

.0001 

1.0  x  10" 3 

S  =  30.4  kx 

s 

-  8  9  .  b  3  7  8 

.0012 

4.58  x  10"5  (NS) 

. 

4 

-89.6379 

.0001 

9.16  x  10"-  (NS) 

Le  par. n't  re  1  nous  perrnet  aussl  d’observer  quo  l'lntensite  au- 
dessus  da  trou  est  on  fait  cello  qu’on  retrouveral  t  en  milieu  homo  gene 
et  quo  les  ray  ms  on  dess^us,  qui  sont  les  rayons  "guides",  possedent 
line  intensite  super  ieure.  I>e  plus,  pour  les  rayons  guides,  l'lntensite 
relative  I  augment*  plus  on  s’eioigne  de  1’antenne.  Ce  dernier  coitpor- 
temont  de  I  nous  procure  un  autre  moyen  de  detecter  le  trou.  Lors  du 
calcul,  en  emettant  les  rayons  avec  un  angle  de  depart  croissant,  on 
pent  reconnoitre  les  rayons  pour  lesquels  l'lntensite  devlent  nettement 
superieure  an  seuil  (i.e  -lO-'  dB)  avant  leur  reflexion  et  identifier 
I'anglo  nil  cela  cesse.  Couch-  il  n'est  plus  requis  d'attendre  que 
1 ' 1 n teg ra 1 1  on  de  terns  les  rayons  soil  complete*  pour  analyser  les 
tables  d'intenslte  resultar.tes  ,  cette  fa<;on  de  proceder  pent  s'averer 
plus  efficace.  Cependant,  er.ant  donne  la  grandeur  des  parametres  en 
je  i  n.ius  devons  et  re  prudents  dans  1 'usage  de  cette  derniere  methode; 
dans  cert. lines  conditions,  les  erreurs  inherentes  au  calcul  ou  a  la 
program.-:  it  ion  pourralent  fausser  les  conclusions. 
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En  connalssant  ainsi  les  rayons  delimitant  le  trou,  il  est  pos¬ 
sible  d'obtenir,  3  l'aide  d'une  technique  iterative  de  lancer  de 
rayons,  la  precision  desirSe  sur  1 ' emplacement  du  trou.  Cependant  on 
salt  qu'en  pratique  un  trou  radio  ne  presente  pas  de  frontieres  si 
abruptes  dues  a  la  diffraction  et  aux  fluctuations  de  l'atmosphere 
qu'on  a  negligees.  Une  precision  de  l'ordre  de  celle  des  fig.  4,  5  et 
6  nous  apparaft  largement  suffisante  pour  les  applications  radars. 

Pour  les  deux  trous  radio  on  peut  done  definir  dans  l'espace  la 
region  "blanche"  ou  l'energie  est  indeterminee  en  optique  geometrique. 
Les  rayons  superieurs  et  inferieurs  d£limitent  le  trou  a  l'origine,  et 
le  second  est  determine  par  son  rayon  superieur  et  la  branche  6upe- 
rieure  de  la  caustique.  II  est  important  qu'un  modele  de  propagation 
utilisant  la  theorie  des  rayons,  et  conqu  pour  evaluer  l'energie  en 
divers  points  de  l'espace,  soit  apte  a  reconnaftre  ces  zones  pour  ne 
pas  chercher  a  atteindre  avec  un  rayon  un  point  situe  a  l'interieur. 

L'intensite  calculee  sur  les  rayons  nous  a  permis  de  constater 
qu'un  conduit  est  di f f ici lenient  decelable  a  partir  de  son  origine; 
l'energie  guidee  a  l'entree  n'etant  pas  superieure  au  cas  normal  de 
faqon  significative.  L'effet  caracteristique  du  conduit  devient  appre¬ 
ciable  a  mesure  que  l'energie  y  est  canalis£e  par  reflexions  multi¬ 
ples. 


6 . 3  Comportement  de  l'intensite  a  la  rencontre  d'une  caustique 

Considerons  encore  une  fois  le  cas  de  la  fig.  6.  Le  rayon  situe 
au  milieu  du  faisceau  reflechi,  et  identifie  ”3"  sur  la  figure,  ren¬ 
contre  la  cuspide  de  front.  La  fig.  7  raontre  l'intensite  calculee  sur 
celui-ci  en  fonction  de  la  distance  parcourue.  Tel  qu'on  pouvait  s'y 
attendre,  on  constate  que  l'intensite  progresse  tres  abruptement  a  la 
cuspide.  Notre  mesure  de  l'intensite  associee  3  un  rayon  apparalt  done 
efficace  pour  deceler  la  presence  d'un  tel  point  singulier.  Pour 
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La  fig.  8  met  en  relief  l'intensite  en  fonction  de  la  position 
dans  le  plan  d-h  pour  le  cas  de  la  fig.  5.  On  y  reconnalt  la  cuspide, 
les  caustiques  superieures  et  inferieures  et  la  zone  d 1  inter ference  due 
aux  trajets  multiples  produite  a  la  reflexion.  Pour  les  points  alimen- 
tes  par  plus  d'un  rayon  une  summation  complexe  a  ete  effectuee  sur  les 
rayons  impliques.  La  longueur  d'onde  choisie  (10  cm)  est  grandement 
inferieure  a  la  resolution  du  graphique.  C'est  pourquoi  les  zones 
d ' interference  apparaissent  brouillees. 


FIGURE  8  -  Intensite  sur  la  surface  S-H  dans  le  cas  du  rayonnement  de 
a  fig.  5 
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7.0  DISCUSSION  ET  CONCLUSION 


Une  theorie  des  rayons  "renouvel^e"  a  §te  recemment  developpee 
au  CRDV.  El le  possede  la  caracteristique  d'integrer,  par  rapport  au 
chemin  optique,  la  trajectoire  des  rayons  et  une  intensite  qui  leur  est 
associee.  Le  present  travail  avait  pour  but  de  verifier  l'applica- 
bilite  de  cette  theorie  a  la  modelisation  de  la  propagation  radar  en 
milieu  heterogene.  II  s'agissait  principalement  de  verifier  dans 
quelle  mesure  1 'information  nouvelle,  apportee  par  le  calcul  d'inten- 
site,  pouvait  contribuer  a  caracteriser  les  phenomenes  "anormaux" 
causes  par  la  refractivite  tels  que  les  trous  radio  et  les  caustiques. 
Les  equations  de  rayon  furent  experimentees  dans  une  situation  de 
conduit  tropospherique  simple  mais  reconnue  comme  representative. 

Cocme  l'effet  de  canalisation  d'energle  est  cause  par  des  rebondisse- 
ments  successifs  de  rayons  rediriges  vers  la  terre  par  la  refractivite 
du  milieu,  il  a  ete  necessaire  de  developper  une  methode  permettant  de 
poursuivre  le  calcul  de  l'intensite  a  la  suite  de  reflexions. 

Un  module  informatique  a  ete  realise  pour  effectuer  l'integra- 
t ion  des  rayons.  Les  resultats  obtenus  ont  ete  exposes  et  analyses  au 
chapitre  6.  On  a  montre  comment,  a  l'aide  du  calcul  de  l'intensite  sur 
un  rayon,  on  peut  facilement  detecter  et  localiser  numeriquement  les 
caustiques  et  cuspides,  et  ce  durant  1 ' in tegrat ion .  Le  reperage  d'un 
trou  radio  produit  3  l'antenne  est  aussi  possible  mais  il  s'est  cepen- 
dant  montre  plus  6ubtil.  En  effet,  a  une  distance  donn£e  de  l'antenne 
lorsqu'on  s'approche  du  trou  d'un  cote  ou  de  l'autre,  on  a  vu  que 
l'intensite  ne  regresse  pas  regul  ierement  de  faqon  significative.  Les 
rayons  se  separent  en  deux  groupes  quasi  homog^nes.  Un  ecart  d'inten- 
site  existe  entre  les  groupes  et  crolt  a  mesure  qu'on  s'eloigne  de 
l'antenne,  mais  il  deraeure  relativemen'  faible  (de  l'ordre  de  10-3  dB); 
c'est  que  la  concentration  d'energie  a  l'entree  du  conduit  avant  le 
premier  rebondissement  est  peu  importante. 
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Puisqu’une  seule  emission  d'un  nombre  restreint  de  rayons  est 
necessaire,  le  calcul  d'lntensite  constitue  un  outil  efficace  pour 
caracteriser  un  patron  de  rayonnement.  Le  reperage  des  zones  singu- 
lieres  peut  permettre  a  un  module  de  propagation  de  decouper  le  milieu 
en  zones  carac te r i s t iques  definies  en  fonction  des  rayons  qui  les 
traversent  et  surtout  de  leur  nombre.  Ce  module  pent  ainsl  evaluer 
rapidement  l'energie  en  un  point  quelconque  sans  tenter  d'atteindre  par 
des  rayons  des  regions  qui  leur  sont  inaccess i bles  (ex:  trous  radio). 

II  connatt  de  merae  les  regions  pour  lesquelles  l'energie  est  indetermi- 
nee  en  optique  geometrique. 

Notons  que  les  observations  et  les  conclusions  qui  ont  ete 
tirees  s'appuienr.  sur  des  resultats  relatifs  plutot  qu'absolus  du 
calcul  de  l'intensite.  La  precision  des  resultats,  que  ce  soit  sur  les 
positions  ou  sur  les  calculs  d'energie,  depend  de  plusieurs  facteurs. 

II  y  a  les  erreurs  inherentes  au  calcul  numerique  qui  dependent  de  la 
precision  du  calculateur,  du  pas  d ' i nt egra t  ion  et  du  nombre  de  rayons 
choisis,  et  aussi  les  erreurs  qui  decoulent  des  approximations  faites  a 
la  reemission.  Bien  entendu  ces  erreurs  d'ordre  numerique  augmentent 
avec  la  distance  franchie.  Une  £tude  pourrait  etre  realisee  pour 
evaluer  la  sens ibi 1 1 te  de  ces  differents  facteurs. 

Les  resultats  obtenus  sont  directement  relies  au  modele  de 
refractivite  adopte  ou,  en  d'autres  mots,  a  notre  connaissance  du 
milieu.  D'une  part,  les  variations  fines  de  l'indice  de  refraction 
(n),  provoquees  principalement  par  la  turbulence,  sont  negligees. 
D'autre  part,  l'indice  de  refraction  ne  nous  est  habituellement  donne 
qu'en  fonction  de  1'elevation  et  on  doit  supposer  qu'il  demeure  cons¬ 
tant  avec  la  distance.  Cette  derniere  hypothese  est  generalement 
valide  sur  une  bonne  distance  mais  au-del3  de  50  km  la  probability  que 
les  conditions  atmospher iques  changent  augmente  considerablement.  De 
plus,  les  caracterist  iques  de  la  surface  ref lechissante  sont  generale¬ 
ment  ma 1  traitees.  La  rugosite  de  la  surface  et,  pa r t i cul ieremen t  en 
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mer,  l'effet  des  vagues  et  leur  mouvement  ne  sont  pas  negligeables .  II 
en  results  que  les  directions  effectives  des  rayons  reflechis  ne  sont 
pas  celles  calculees.  Encore  une  fois,  les  erreurs  engendr£es  sont 
cur/.: la 1 1  ves  de  rebondissement  en  rebondissement .  Toutes  ces  impreci- 
sions  sur  l'indice  de  rlfraction  et  la  surface  contribuent  a  brouiller 
le  schema  statique  calcule.  Ceci  explique  pourquoi,  en  pratique,  meme 
pour  A  petit,  l'energie  mesuree  dans  un  trou  radio  est  non  nulle  et  les 
caustiques  ne  sont  pas  si  clairement  definies. 

Le  nombre  de  rayons  a  emettre  et  le  nombre  de  zones  caracteris- 
tiques  3  considerer  pour  evaluer  l'energie  augmentent  rapidement  3 
mesure  que  l'on  s'£loigne  de  l'antenne.  Aussi,  comme  les  erreurs 
decrites  ci-haut  croissent  avec  la  distance,  la  theorie  des  rayons 
devient  moins  efficace  a  de  grandes  distances  de  l'antenne.  II  peut 
etre  preferable  dans  ce  cas  de  faire  appel  a  la  theorie  des  modes,  meme 
si  cela  implique  1 'utilisation  de  formules  empiriques.  L'usage  de  la 
theorie  des  modes  est  en  effet  plus  pertinent  3  des  distances  ou  le 
nombre  de  modes  toujours  persistants  est  2  ou  3  au  maximum  (ggnerale- 
ment  3  partir  d'environ  100  km).  Par  contre  sur  de  plus  courtes  dis¬ 
tances  le  nombre  de  rebondissements  a  traiter  est  reduit  et  l'approche 
des  rayons  devient  beaucoup  plus  pratique  vu  la  complexite  croissante 
de  la  theorie  des  modes.  II  demeure  toutefois  qu'au  voisinage  des 
caustiques,  comme  on  l'a  mentionne  au  chapitre  3,  la  theorie  des  rayons 
ne  tient  plus  et  on  doit  recourir  a  l'€lectromagnetisme  classique  si  on 
desire  resoudre  plus  finement  le  probl3me  de  l'energie. 

Enfin,  pour  etre  plus  representatif  et  plus  juste,  un  module  de 
propagation  devrait  tenir  compte  des  autres  facteurs  qui  affectent 
consi derablement  la  propagation  des  ondes  61ectromagn£tiques .  En  effet, 
il  n'est  pas  realiste  de  traiter  isolement  les  effets  de  la  refracti- 
vlte  sans  tenir  compte  aussi  des  autres  phenomenes  causes  par  les  con¬ 
ditions  atraospher iques  telles  l'absorption  moleculaire  et  1 ' at tenuat ion 
dues  aux  hydrome teores .  On  doit  aussi  pouvoir  faire  intervenir  un  dia- 
gramme  de  rayonneraent  d'antenne  reel. 
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6.0  TRAVAI X  FlTL'KS 

Pours.iivant  noire  object  if  de  paranet  riser  la  propagation  en 
milieu  heterogene,  nous  projet  ms  d'etudier  le  cas  de  profils  de 
refractivite  modelises  a  partir  de  donnees  me  teorologiques  et  utilisant 
des  criteres  de  continuity  sur  la  fonction  N  et  ses  premiere  et  seconde 
derivees  (ref.  23).  Nous  conptons  analyser  le  comportement  de  notre 
nodele  de  propagation  dans  ces  conditions  dans  le  but  de  generaliser 
les  conclusions  tirees  des  resultats  obtenus  avec  les  profils  de 
Boithias.  Cela  nous  permettra  aussi  de  verifier  les  consequences  sur 
le  rayonneaent  des  hypotheses  posees  lors  de  la  modelisation  de  la 
fonction  de  refractivite.  Nous  pourrons  ainsi  mettre  au  point  une 
modelisation  adequate  des  profils. 
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ANNEXE  A 


Equation  eikonal  d£rivee  des  Equations  de  Maxwell 
dans  le  cas  ou  la  polarisation  est  lineaire  et  invariante 

Dans  un  milieu  sans  charge  ni  courant,  les  equations  de  Maxwell 


V«e  E 

=  0 

[A.!] 

V.yH 

=  0 

[  A  -  2  ] 

-1  3yH  | 

E  = 

c  3 1 

.  A.  3j 

1  3cE  | 

H  = 

A .  4  J 

c  3t 

■f  4 

ou  E  et  H  sont  les  fonctions  vectorielles  de  champs  electrique  et 
magnetique.  c  est  une  constante  egale  a  la  vitesse  de  la  lumiere  dans 
le  vide,  et  e  la  constante  dielectrique  et  y  la  permittivite  sont  des 
parametres  carac t£risant  entierement  le  milieu  de  propagation.  Ces 
derniers  sont  reels  dans  le  cas  ou  le  milieu  est  sans  perte  (milieu 
isotropique) .  Ils  sont  relics  a  l’indice  de  refraction  n  selon: 


n  * 


A. 5] 


Considerons  l'expresslon  du  champ  monochroma t ique 
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les  equations  A.l  a  A. 4  deviennent: 


V*  £  E 

=  0 

00 

cu 

0 

r  T 

V*  uH 

=  0 

A. 8b 

O  L 

+ 

[ A. 8c  ] 

V  x  E  = 
0 

jk0uH 

-*■ 

-► 

r  q 

V  x  H  = 
0 

-  jk  cH 

J  0 

(A.8d 1 

Si  la  variation  temporelle  de  £  et  u  est  negligee  au  cours  de  la 

*► 

propagation  (considerant  un  milieu  statique),  en  substituant  dans 
[a. 8a:  on  trouve: 

E^* V$  =  0  [A. 9] 


E  *7e 

et  V*E  =  -  — r  ■—  =  -E  »Vlne 


De  la  meme  faqon,  on  trouve  a  l'aide  de  [ A  -  8b ] : 


H  •  Vq>  =  0 

r 


et  V*H  =  H  *Vlnr 
r  r 


En  6ubstituant  maintenant  E  dans  [a.8c  ,  on  aura: 


V  x  E  +  jk  (V4,  x  E  )  =  jk  pH 


I  a. 10] 


a.h  ] 


a. 12 


i  A. 13] 
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Le  premier  terme  etant  nul  (equation  A. 9),  on  obtient  l'equation 
eikonal: 


2  =  n2 


A. 21 


II  est  aussi  possible  d'obtenir  le  meme  resultat  a  partir  de 
l'equation  d'onde  en  milieu  isotropique  et  homogene,  soit: 


v2e  -  nL  lii  =  _IiL  x  (V  x  e)  -  v.  (li-E) 
T7  at2  w  e 


A. 22 


sans  faire  1 'approximation  Vp  =  0.  Le  calcul  n'est  pas  pr£sente  ici 
II  est  beaucoup  plus  long  et  comae  il  ne  s'agissait  que  d'exposer  le 
principe,  on  a  juge  preferable  de  deriver  l'equation  eikonal  par  une 
methode  plus  directe. 


^  m-'r% 


A 


T0-A169  553  HETHODE  OES  RAVONS  AVEC  CALCUL  D'  INTEN5ITE  flPPLIOUEE  A  M 
LA  PROPAGATION  ANO  (U)  DEFENCE  RESEARCH  ESTABLISHMENT 
VAL CARTIER  (OUEBEC)  D  DION  JUN  86  DREV-4487/86 
UNCLASSIFIED  F/G  28/14  NL 


DREV 

R-4407/86 

UNCLASSIFIED 

CRDV  R-4407/86 

FILE: 

3632D-005 

DOSSIER:  3632D-005 

JUNE 

1986 

JUIN  1986 

METHODE  DES  RAYONS  AVEC  CALCUL  D’INTENSITE 
APPLIQUEE  A  LA  PROPAGATION  "ANORMALE" 


par 

D.  Dion 


CORRIGENDUM  -  ERRATA 


Les  iquations  et  le  texte  ci-joints  remplacent  ceux  qui 
apparalssent  dans  le  CRDV  R-4407/86.  Ils  sont  Imprimis  sur  du 
papier  qui  devlent  adhisif  lorsqu'il  est  mouill§. 


The  attached  equations  and  text  supersede  those  appearing  In  DREV 
R-4407/86.  They  have  been  printed  on  paper  that  is  adhesive  when 
wet. 


Quebec,  Canada 


November /novembre  1986 


SANS  CLASSIFICATION 


ERRATA 


i 

i 

Page  3:  Liquation  [  2 . 1 J  doit  i tre  remplac^e  par: 


(70(x,  y,  z))2  -  n2(x,  y,  z)  [2.l] 


et  l’equation  [2.2]  par: 


-► 


Page  5:  Le  troisieme  paragraphe  doit  etre  remplace  par: 


On  montre  dans  la  r€f .  4  que  l'intensite  le  long  d'un  rayon 
obeit  £  la  relation: 


!(♦)  a 


1 

- T - 

n(a  x 


b)-p 


[2.8] 


fcn  pratique,  lorsque  les  rayons  se  croisent,  a  la  suite  d'une 
reflexion  ou  en  traversant  une  caustique,  a  ou  b  changeront  de'direc- 
tion,  done  de  signe.  On  assiste  en  fait  I  un  renversement  de  l’onde 
sur  elle-m?me.  Bien  entendu,  afin  d’Sviter  d'obtenir  des  intensites 
negatives,  il  est  necessaire  de  considerer  dans  les  calculs  la  valeur 
absolue  de  l'expression  [ 2.tt]  . 


Page  10:  On  doit  remplacer  le  premier  paragraphe  par: 


Le  gain  en  amplitude  resultant  variera  ainsi  entre  0  et  2  A.  ;  ce  qui 
correspond  a  un  gain  en  intensite  allant  de  0  a  (I  A^)2.  Lorsqu'un 
rayon  traverse  une  caustique,  on  s'aperi;oit  qu'il  est  tourne  sur 
lui-meme.  Ceci  est  indique  dans  not  re  modele  par  le  changement 
d  orientation  des  vecteurs  a  ou  b.  On  assisterait  done  2  un  dephasage 
de  jt/2.  Les  references  2  et  5  predisent  en  effet  ce  dephasage.  II  est 
a  remarquer  que  ce  phenomeme  est  en  quelque  sorte  aussi  observe  lors  de 
la  reflexion  d'une  onde  Slectromagnetique  a  polarisation  verticale. 


Page  20: 


A  la  fin  de  la  premiere  phrase,  on  doit  lire  au  lieu  de  "... 
longueur  optique  0^. 


longueur  optique  4^. 


Page  24:  Sur  le  graphique  de  la  figure  3,  on  doit: 
1)  remplacer  -0,04  N/km  par: 


-40  N/km 


2)  remplacer  -0,5  N/km  par: 


-500  N/km 


rn 

a 


Page  42:  Au  milieu  de  la  page,  on  doit  reiaplacer  "En  substituant 
(A. 14]  dans  [  A  - 15]  ,  on  aura:" 


En  substituant  [A.14]  dans  [a.16],  on  aura: 


Page  43:  On  doit  remplacer  le  debut  du  dernier  paragraphe  par: 


11  est  aussi  possible  d'obtenir  le  m§me  resultat  a  partir  de 
l'equation  d'onde  en  milieu  isotropique,  soit: 


